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We investigate tiieoretically the effects of intrinsic spin-relaxation on the spin-dependent transport 
through a single-molecule magnet (SMM), which is weakly coupled to ferromagnetic leads. The 
tunnel magnetoresistance (TMR) is obtained by means of the rate-equation approach including 
not only the sequential but also the cotunneling processes. It is shown that the TMR is strongly 
suppressed by the fast spin-relaxation in the sequential region and can vary from a large positive to 
slight negative value in the cotunneling region. Moreover, with an external magnetic field along the 
easy-axis of SMM, a large negative TMR is found when the relaxation strength increases. Finally, 
in the high bias voltage limit the TMR for the negative bias is slightly larger than its characteristic 
value of the sequential region, however it can become negative for the positive bias caused by the 
fast spin-relaxation. 

PACS numbers: 75.50.Xx, 75.76.-fj, 85.65.-f-h, 85.75.d 



I. INTRODUCTION 



The electron transport through magnetic molecules 
has attracted much attention recently in molecular 
spintronicsi"— , because of its potential applications in 
information storage and processing. The single-molecule 
magnet (SMM), which has a high intrinsic spin and an 
easy-axis magnetic anisotropy, is of particular interest. 
Some peculiar phenomena such as Coulomb-blockade, 
magnetic excitations, complete current-suppression, and 
negative differential conductance have been observed ex- 
perimentally in the quantum transport through a SMIV^. 
The Berry phase'^ and Kondo effects^ on the average 
current (or differential conductance) are predicted the- 
oretically. Moreover, the full counting statistics (current 
fluctuation)^ and spin fluctuation of magnetic molecules^ 
are also studied. On the other hand, the magnetic 
switching of SMM is realized by thermal spin-transfer 
torque^, spin-bias^, and spin-polarized current injected 
from ferromagnetic (FM) electrodes^iii"— . It is pre- 
dicted that the spin-current polarization can be reversed 
through a SMM with FM leads^^ and even a pure spin- 
current can be generated in normal leads by the thermo- 
electric effectsi^. Spin-diode behavior— and memristive 
properties^? are also observed in spin-polarized trans- 
port through a SMM. The current difference between the 
parallel (P) and antiparallel (AP) configurations of two 
FM electrodes known as the tunnel magnetoresistance 
(TMR)2*. was systematically investigated for the SMM by 
Misiorny et al. in the sequential, cotunneling, and Kondo 
regionsiSr— . The TMR is positive usually, since the P 
current is greater than the AP one, while a large negative 
TMR is shown to exist in the SMMiS. Furthermore, the 
tunneling anisotropic magnetoresistance (TAMR) with 
tunable magnitude and sign is found in SMM junctions 
with only one FM electrode^^. 

A long spin-relaxation time of the SMM is required in 



the quantum information processings^ and information 
storagei^,. However, in addition to the spin-relaxation 
induced by the electron tunneling through molecules, 
the SMM itself may have intrinsic spin-relaxatioE^i2i2^, 
which originates from the interactions with environmen- 
tal spins. A short spin- relaxation time is found in a recent 
experiments^ in which a spin-polarized scanning tunnel- 
ing microscopy is used to pump the spins of Mn atoms. 
Moreover, the experimentally observed quenching of con- 
ductance steps2^ has been successfully explained in terms 
of the spin-relaxation2£. The effects of intrinsic spin- 
relaxation on the current-induced magnetic switching of 
SMM have been studied theoreticallyi^. When the relax- 
ation of SMM is very fast, the fine structure of transport 
is suppressed, but a high spin-polarized current in non- 
magnetic leads can be generated with an external mag- 
netic field^^. For the quantum-dot (QD) systems, on the 
other hand, the intrinsic spin-relaxation can result in the 
TMR suppression'^ - and reversal^ in the sequential and 
cotunneling regions, respectively. It is also shown that 
the peak heights of differential conductance in QDs are 
increased by the spin-relaxation5^. 



In this paper we study the spin-polarized transport 
through a SMM with FM leads in both sequential and 
cotunneling regions and pay particular attention on the 
intrinsic spin-relaxation, which affects remarkably the 
transport properties. A large negative TMR can be ob- 
tained when a magnetic field is applied along the easy- 
axis of SMM. Sec. II is devoted to the rate-equation for- 
malism for the spin-polarized transport through a SMM 
and the related TMR. We demonstrate and discuss the 
effects of intrinsic spin-relaxation on the current, differ- 
ential conductance, magnetization and TMR based on 
the numerical results in Sec. III. Finally, the conclusion is 
given in Sec. IV. 
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FIG. 1: (Color online) Schematic diagram of electron trans- 
port through a SMM weakly coupled to ferromagnetic elec- 
trodes. The magnetizations of the two leads are coUinear with 
the magnetic easy axis of SMM (as z axis). The left (L) and 
right (R) electrodes are connected with the bias voltage V/2 
and -V/2, respectively. 



II. MODEL AND METHOD 

The schematic diagram of electron transport through 
a SMM coupled to two external FM leads is shown in 
Fig. 1. It is assumed that the lead- magnetizations are 
collinear with the magnetic easy-axis of SMM in either 
P or AP configuration. The total Hamiltonian of the 
system can be written a a^^i^^i^^i^"^ 



H 



H, 



SMM 



Hj 



(1) 



The first term Hieads = Z]a=L,Rl]kcr ^akcJ^ktrCako- de- 
scribes noninteracting electrons in electrodes, where cj^j.^^ 
(cakcr) is the creation (annihilation) operator for an elec- 
tron with wave vector k and spin a in the lead a, and 
£aka is the corresponding electron energy. The spin 
polarization of FM lead a is defined as Pa = {pa+ — 
Pa-)/{pa+ + Pa-), with Pa+{-) denoting the density of 
states for the majority (minority) electrons. The second 
term 



HsMM — /_^ 



= ^(e - eVg)dld„ + Ud\d-^dld^ - Js • S 



K{S'-f - B{3'' + S'') 



(2) 



denotes the SMM Hamiltonian, where the operator 
d^ (do-) creates (annihilates) an electron in the low- 
est unoccupied molecular orbital (LUMO) level and s= 
Xlcrcr' d\{(Jaa' l'^)da' is the electron spin operator with a 
being the vector of Pauli matrices. The single-electron 
energy e in the LUMO level is tunable by gate voltage 
Vg . The Coulomb energy denoted by U is for the double 
occupancy of the LUMO level. The exchange coupling 
J between the local spin of molecule S and the LUMO 
electron-spin s can be either FM (J > 0) or antiferromag- 
netic (AFM) (J < 0). The easy-axis anisotropy of the 
SMM is characterized by the parameter K [K > 0) and 
an external magnetic field B is applied along the easy- 
axis, where B contains the factor qijlb- Since the SMM 
Hamiltonian Hsmm commutes with z component S^ of 
the total spin operator St= s -I- S, it can be easily di- 
agonalized in terms of molecular many-body eigenstates 



|n, S't;m) with n denoting the charge state, St the total 
spin quantum- number, and m the eigenvalues of S^ . 

The tunneling Hamiltonian between the LUMO level 
and leads isHiiSiS^iSi 
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where the parameter t^ is the tunneling coupling- 
constant between the molecule and lead a, and the corre- 
sponding spin-dependent tunnel-coupling strength is de- 
noted by Taa = 27r/5Q,cr \ta\ • Furthermore, in terms of 
the spin polarization p„ of lead a, we can rewrite the 
tunnel-coupling strength as Ta± = ^a{^ ± Va)l^ for 
the spin-majority (spin-minority) electrons of the lead 
a, with Va = ra+ + Fq-. 

Since the SMM-electrode couphng is assumed to be 
sufficiently weak, i.e., Fcto- <^ ksT, the molecule relaxes 
rapidly to the eigenstates of Hsmm (rapid dephasing)^^. 
Therefore, the rate-equation approach can be adopted 
to study spin-dependent transport through the SMM. 
Here, we consider phenomenologically the intrinsic spin- 
relaxation processes, which may drive the SMM to an 
eigenstate of lower energy. The transport dynamics is 
well described by the following rate equation. 
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with Pi denoting the population probability of the SMM 
state \i). The sequential rate W" describes the transi- 
tion from the molecule state \i) to \i') along with the elec- 
tron tunneling into or out of lead a, and W^^^, denotes 
the cotunneling rate from lead a to a'. The above sequen- 
tial and cotunneling rates can be calculated in terms of 
the T-matrix with the help of generalized Fermi golden- 
rulei^i^"— . Moreover, the transition selection-rules of 
spin relaxation are An = and Am = ±1, which are sim- 
ilar to the cotunneling processes. The phenomenological 
spin-relaxation rate from the state \i) to \i') readsi^i^ 
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where Trei is the relaxation time and Si denotes the en- 
ergy of state \i). In addition, a dimensionless parameter 
7 = h/ {TTrei) called the spin- relaxation strength is intro- 
duced to compare Trei with the sequential tunneling time 
h/T. 

For the steady-state transport we have the stationary 
probabilities with the condition ^^ = 0. Therefore, the 
stationary current through the lead a is given by 
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where the current flows from the left electrode to right 
one, and the magnetization of SMM is (S^) = '^■rriiPi. 
The TMR is defined as 



TMR 



TMR: 



Ip — Iap 
Iap 



(7) 



where Ip {Iap) denotes the current through the SMM in 
the P (AP) configuration. 



III. RESULTS AND DISCUSSION 

We now analyze the effects of intrinsic spin-relaxation 
in both the sequential and cotunneling regions with the 
assumption that the bias voltage is symmetric at the 
SMM-electrode tunnel junction, i.e., fiL — —fJ-p — V/2. 
In the following numerical calculations, the parameters 
are chosen as 5 = 2, e = 0.5 meV, \J\ ~ 0.2 meV, 
U = 1 meV, K = 0.05 meV, F = F^ = F^ =:= 0.001 
meV, PL =Pr = 0.5 and fc^T = 0.04 meV. 

The TMR value indicated by color index as a func- 
tion of the gate {Vg ) and bias {V) voltages is plotted 
in Fig. 2. The electron transports are dominated by co- 
tunneling transitions in the voltage-regions labeled with 
the numbers n = 0, 1, 2. While in the rest of regions 
the transports are dominated by sequential processes. 
The left panels [(a)-(c)] of Fig. 2 are referred to the 
results without the intrinsic spin-relaxation as a com- 
parison, and the corresponding TMR for the relaxation 
strength 7 = 1 is shown on the right panels [(d)-(f)]. 
Comparing Figs. 2(a)-2(c) and Figs. 2(d)-2(f), we can 
see the obvious variation of TMR induced by the intrin- 
sic spin-relaxation, which in most cases leads to the TMR 
suppression. However, when an external magnetic field 
is applied along the easy-axis of SMM, the TMR can be 
increased seen from the Figs. 2(c) and (f). In particu- 
lar, a negative TMR value extends from cotunneling to 
sequential regions. 

The bias voltage dependence of current, differential 
conductance, magnetization of SMM, and TMR are 
shown in Fig. 3 for different relaxation strength 7 
with the gate voltage Vg = —0.3 mV, at which the 
ground state of SMM is |0, 2;±2). For the P con- 
figuration of FM leads, the conductance peak-1 [Fig. 
3(b)], which corresponds to the transition |0, 2;±2) <^ 
|1, 5/2; ±5/2), appears with the increase of bias voltage 
V and is not affected by the intrinsic spin-relaxation 
since the transition is between the ground states of 
n = and n = 1. When the initial state is excited, 
the corresponding peak-2 (|0,2;±1) <=> |1, 3/2; ±3/2)) 
and peak-4 (|1, 3/2; ±3/2) ^ |2,2;±2)) disappear?-? in 
the existence of spin-relaxation. However, the peak-3 
(|0, 2; ±2) O |1, 3/2; ±3/2)) still exists and especially the 
height of peak-5 (|1, 5/2; ±5/2) ^ |2,2;±2)) increases 
along with the population probability of the ground state 
|1, 5/2; ±5/2). This enhancement phenomenon results 
from the relaxation induced decay of excited states. In 
consequence, the current Ip between the peak-1 and 
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FIG. 2: (Color online) TMR as a function of the bias and gate 
voltages in the case of FM coupling [(a),(d)], AFM coupling 
[(b), (e)], and FM coupling with a magnetic field along the 
easy-axis of SMM [(c) , (f )] for different relaxation strength 7 — 
(left panel) and 7 = 1 (right panel). The parameters are: 
S = 2, e = 0.5 meV, |J| = 0.2 meV, U = 1 meV, K = 0.05 
meV, ksT = 0.04 meV, pL = pn = 0.5, F = Fl = Fa = 0.001 
meV, and Jo — 2eT /h. 



peak-3 is slightly increased, while the current between 
the peak-3 and peak-5 is decreased as shown in the Fig. 
3(a). When the bias voltage is high enough such that 
all channels enter into the transport window, the cur- 
rent Ip remains in a constant value independent of the 
relaxation. In addition, the magnetization of SMM van- 
ishes when the relaxation strength increases to 7 = lO'^ 
as shown in Fig. 3(c) (dotted line). The transport prop- 
erties of AP configuration are quite different from the P 
case. For example, in the absence of spin-relaxation the 
magnetization of SMM [Fig. 3(g)] tends to a large nega- 
tive value resulting from the electron spin-flip processi^, 
and the steady transport is dominated by the negative- 
eigenvalue states oi S^ . The spin- relaxation, however, re- 
moves the above tendency and more positive-eigenvalue 
states contribute to the transport. Therefore, the current 
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FIG. 3: (Color online) The bias voltage dependence of current 
I, differential conductance G, and the magnetization of the 
SMM for the parallel [left panel (a)-(c)] and antiparallel [right 
panel (e)-(g)] configurations with the gate voltage Vg — —0.3 
mV. TMR plots are shown in (d) and (h) for the FM and 
AFM exchange couplings, respectively. 
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FIG. 4: (Color online) The bias voltage dependence of Ip (a), 
Iap (e), GP (b), Gap (f), {S!)p (c) {S!)^p (g), and TMR 
for the FM [(d)], AFM [(h)] exchange couplings with the gate 
voltage Vg = 1 mV. The cotunneling currents Ip and Iap are 
shown in the insets of (a) and (e), respectively. 



Iap [Fig. 3(e)] and differential conductance Gap [ Fig. 
3(f)] increase monotonously with the relaxation strength 
7. The 7-dependence of TMR spectrum for both cases of 
FM and AFM exchange couplings are shown in Fig. 3 (d) 
and (h), respectively. At the low bias voltage, electron 
transport is dominated by elastic cotunneling processes 
via ground states of n = and the TMR is essentially not 
affected by the relaxation. With the increase of bias volt- 
age, the inelastic cotunneling plays a role in the transport 
and the current Iap increases due to the intrinsic spin- 
relaxation. Therefore, the TMR almost vanishes at the 
large value of 7. In the sequential tunneling region, the 
fast relaxation also leads to vanishing TMR. In addition, 
the energy of state |1, 3/2; to) for the AFM exchange cou- 
pling is lower than that of the state |1, 5/2; m), and the 
TMR decreases (with the increase of relaxation strength) 
faster compared with the FM case. 

In Figs. 4(a)-4(c), we display the current, differential 
conductance and magnetization of SMM in the P config- 
uration as a function of the bias voltage V for the FM 
exchange coupling and Vg — 1 mV. In this case, the gate 
voltage corresponds to the electron- hole symmetry point, 
thus the magnetization of SMM is always zero indepen- 



dent of the bias voltage V [Fig. 4(c)] and the spectrum 
of differential conductance has only one peak seen from 
Fig. 4(b). Moreover, the peak-position has a slight shift 
in the existence of intrinsic spin-relaxation due to the in- 
crease of population probability of state ]1, 5/2; ±5/2). 
For the P configuration, it can be seen from Fig. 4(a) 
that Ip is suppressed by the spin-relaxation in both the 
cotunneling [the inset of Fig. 4(a)] and sequential tun- 
neling regions. But for the AP configuration, the current 
[Fig. 4(e)] and differential conductance [Fig. 4(f)] are 
enhanced by the relaxation. In the Fig. 4 (d), we show 
the bias voltage dependence of TMR with FM exchange 
coupling for different relaxation strengths. The intrinsic 
spin- relaxation leads to the decrease of TMR from a large 
positive to small negative value in the cotunneling region. 
This can be understood from the behaviors of Ip and Iap 
[see the insets of Figs. 4(a) and 4(e)]. For the AFM ex- 
change coupling, the TMR exhibits a different behavior 
from the FM case. At low bias voltage, in which elec- 
tron transport is dominated by cotunneling processes via 
ground states |1, 3/2; ±3/2), the intrinsic spin- relaxation 
has a little effect on the TMR. When the bias voltage 
increases, the enhancement of TMR induced by the spin- 
relaxation becomes obvious. However, the large negative 



TMR almost disappears due to the fast relaxation [Fig. 
4(h)]. The above behavior of TMR for the AFM coupling 
is due to the decrease of Iap in the cotunneling region. 
When all of transport channels are open in the sequen- 
tial region, the TMR of both FM and AFM exchange 
couplings approaches zero when the relaxation strength 
increases to a sufficiently large value. 

Figure 5 shows the effects of the intrinsic spin- 
relaxation on the bias voltage dependence of TMR in the 
existence of a magnetic field, which is applied along the 
easy-axis of SMM. The longitudinal magnetic field lifts 
the degeneracy of states |n, 5't;±m) and the fast spin- 
relaxation leads to the transition from the state |n, St', m) 
to the ground state \n, St', St) during the transport pro- 
cess. Moreover, this magnetic field also leads to the 
symmetry-breaking with respect to the bias reversal in 
the AP configuration, but the current remains symmet- 
ric in the P configuration. The TMR is illustrated in 
Fig. 5(a) for V^ = 2.3 mV, where the transport is mainly 
through the ground state |2,2;2) in the cotunneling re- 
gion. It is found that the spin-relaxation induced TMR- 
variation is negligibly small in the negative bias voltage, 
while the spin-relaxation can lead to a great change of 
TMR varying from positive to negative values in the pos- 
itive bias voltage. Since the currents Ip for both posi- 
tive and negative bias voltages are suppressed equally 
by the spin-relaxation in the P configuration, the above 
asymmetric behavior of TMR is determined by the trans- 
port in the AP configuration only. For this case, there 
exists a competition between the spin-flip process^^ (in- 
duced by the inelastic cotunneling) and the intrinsic spin- 
relaxation. At the negative bias voltage, the above two 
mechanisms make the SMM be trapped in the ground 
state |2,2;2) and the current is constant. However, at 
the positive bias voltage, the spin-flip process makes the 
SMM tend to the state |2, 2; — 2), but the spin- relaxation 
drives it to the state |2,2;2). In consequence, the ef- 
fect of relaxation can compensate the decrease of Iap 
induced by the spin-flip. Moreover, the TMR for Vg = 1 
mV [Fig. 5(b)] and Vg = -0.3 mV [Fig. 5(c)] exhibits a 
similar behavior in the cotunneling region. 

In order to understand the TMR behavior in the se- 
quential region at V^ = 2.3 mV [see Fig. 5(a)], Ip and 
Iap are plotted in Fig. 5(d) for the relaxation strength 
7 = and 7 = 10'^. At the negative bias voltage, both 
the spin-flip processiiii^ii^'i^i^ and the intrinsic spin- 
relaxation can increase the population probabilities of 
large positive-eigenvalue states of S^ so that the relax- 
ation induced change of Iap is very small. However, at 
the positive bias voltage, the current Iap is increased 
much by the relaxation. On the other hand, for both 
positive and negative bias voltages in the P configura- 
tion, the spin-relaxation results in the enhancement of 
sequential transition |2, 2; 2) <^ |1, 5/2; 5/2) for the spin- 
down electrons, which are minority in FM electrodes, and 
thus Ip decreases. As a result, the corresponding TMR 
can vary from positive to negative values along with the 
increase of the relaxation strength at a fixed bias volt- 




FIG. 5: (Color online) TMR as a function of bias voltage 
for different gate voltages (a) Vg = 2.3 meV (b)l meV (c)0.3 
meV. (d) Plots of absolute currents for different relaxation 
strengths with a magnetic field along the easy-axis of SMM. 



age. In particular, a large negative TMR is generated 
by the fast spin-relaxation at the positive bias voltage. 
Fig. 5(c) shows the TMR as a function of the bias volt- 
age for Vg = —0.3 mV, where the ground state of SMM 
is [0,2; 2). Different from the Vg — 2.3 mV case, the 
TMR is enhanced by the spin-relaxation at the nega- 
tive bias voltage while it is suppressed even to zero at 
the positive bias voltage. This is because that the in- 
trinsic spin-relaxation leads to the increase of Ip by the 
enhancement of majority-electron transport through the 
sequential transition |0, 2; 2) <^ |1, 5/2; 5/2) in the P con- 
figuration. In addition, when the bias voltage is high 
enough such that all transport channels are open, the 
Iap is slightly reduced by the spin-relaxation at the neg- 
ative bias, but it can be greatly enhanced at the positive 
bias. Therefore, the TMR is slightly above the charac- 
teristic value of the sequential region in the negative bias 
and becomes negative in the positive bias seen from the 
Figs. 5(a)-5(c). 



IV. CONCLUSION 

In summary, the impact of the intrinsic spin-relaxation 
on electron transport through a SMM with FM leads is 
investigated by the explicit calculation of spin-relaxation 
dependence of transport quantities such as current, dif- 



ferential conductance, magnetization and TMR in both 
the sequential and cotunneling regions. The differen- 
tial conductance peaks can be completely suppressed by 
the spin-relaxation for the transition from initial excited- 
states in the P configuration. The magnetization varies 
from a large negative value to zero with the increase of 
spin-relaxation strength in the AP configuration, and the 
current Iap is enhanced. Therefore, the TMR is com- 
pletely suppressed in the sequential region and can vary 
from a large positive to slight negative value in the co- 
tunneling region. When an external magnetic field is ap- 
plied along the easy-axis of SMM, the TMR is asymmet- 
ric with respect to the bias reversal. For the cotunneling 
process, the TMR is much more sensitive to the variation 
of spin-relaxation in positive bias than in negative bias, 
especially, it can change from positive to negative values 



in the former case. In the sequential region, a large nega- 
tive TMR value can be generated by the spin-relaxation 
in positive bias. Finally, in the high bias limit the TMR 
at negative bias is slightly larger than the typical sequen- 
tial TMR value, while at positive bias it becomes negative 
due to the spin-relaxation. It is expected that the above 
transport properties can be observed experimentally in 
FM-SMM-FM spin valves^. 
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